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Abstract 



A reanalysis of identified muon neutrino interactions from 1MB 3 yields 
bounds on sin (20) and Am 2 . The limit sin 2 (20) < 0.72 is in conflict with 
the recent announcement of a neutrino mass. 

Subject headings: Cosmic Rays — Elementary Particles — Neutrino Os- 
cillations 
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Atmospheric neutrinos have always held the possibility of providing a large range of 
propagation coupled with a near source to provide a convenient calibration. With a reason- 
able choice of angular cuts a sample of neutrinos traveling about 10,000 km can easily be 
compared with a sample traveling on the order of 20 km collected in the same experiment 
at the same time. Early work ]I] using these concepts was hampered by low statistics and 
poor particle identification. In [|lj] a mixed sample of 25 upward neutrino interactions was 
compared with a sample of 25 downward neutrino interactions. The neutrino propagation 
distance is a function of the neutrino direction, so scattering can, in principle, make the 
distance estimate ambiguous. In reality the propagation distance varies very slowly with 
angle near the vertical so as long as regions near the vertical are compared scattering will 
have a negligible effect on the sensitivity. In most analyses regions containing 20% of the 
solid angle over the upward and downward direction have been used. 

Depending on the location of the detector the upward and downward fluxes may be 
influenced differently by geomagnetic effects . 1MB was located at 52° north geomagnetic 
latitude. This meant that the local flux, coming from above, had a lower geomagnetic cut off 
than the Earth at large. As such one would expect a greater rate of downward going events 
since the lower cut off permitted more of the extraterrestrial cosmic ray flux to descend and 
hit the atmosphere. But such effects are small and have not been observed in a statistically 
significant way even in a sample of 401 neutrino interactions 0. In the 401 1MB 1 sample 
55 events were measured in the upward going 20% of the solid angle and 65 were found in 
a comparable portion of downward solid angle, comparison of the 15 upward events and 21 
downward events which were classified as muon neutrino induced because of the presence of 
a muon decay in the event yielded limits on Am 2 in the range of 5 x 10~ 5 eV 2 . 

1MB 3 provided a comparable sample of events but with better light collection and with 
morphologically based muon identification methods || as well as those based on muon decay 
signatures. An analysis comparing the shape of the energy distribution for the upward and 
downward 20% of the solid angle M confirmed the 1MB 1 results but with a larger sample 
and with particle identification. The 1MB 3 work was done with 34 downward going and 
32 upward going v^. These two samples have energy distributions which are statistically 
indist inguishable . 

The excluded region derived in || is bounded in the range 4.5 x 10~ 5 < Am 2 < 1.1 x 10~ 4 
eV 2 . Outside of this range distortion of the upward going neutrino spectrum by neutrino os- 
cillations would not make its shape statistically distinguishable from the measured downward 
spectrum. 

The contained neutrino interactions can be used to study neutrino oscillations over a 
broader range of Am 2 . The lowest limits on Am 2 are found in comparisons of the shape of 
the observed contained event energy spectra for samples of events near the vertical. Such 
a test is insensitive to detailed flux calculations, since the downward sample constitutes a 
measured near source. The method has a limited range because for larger Am 2 the upward 
going sample would be fully mixed via multiple oscillations and hence would have a spectrum 
shape that was similar to the unoscillated downward sample. The upward sample would be 
reduced by an overall factor, which is 1 — |sin 2 (2#). So comparison of rate of upward to 
downward going interactions can extend the sensitive range. This comparison eventually 
fails at even higher Am 2 because at large Am 2 both the upward and downward going samples 
will have oscillated. In such cases oscillations may be noted by deviations of the observations 
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from expectations. Estimates of Am 2 would be unreliable since no distance scale would be 
present in the observations. 

The first two methods, comparison of up and down spectral shape and comparison of 
up and down interaction rate are insensitive to most systematic errors in the neutrino flux 
calculations. For example they are insensitive to the normalization uncertainty, the to 
v e ratio uncertainty and the temporal modulation caused by the solar wind's impact on 
the Earth's magnetosphere. They may be influenced by geomagnetic effects that produce 
modest variations of the neutrino flux at various points on Earth. The downward neutrino 
component is produced locally and so is sensitive to local magnetic field properties. 

The result of can be extended by comparing the rate of upward going to downward 
going interactions. For the range of Am 2 of greatest interest for the Super Kamioka f| 
results one would expect a two to one rate difference 0. 

For the data of reference 0] the up to down event ratio is 
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— = 0.94 ±0.23 > 0.64 
34 

where the bound is the 90% confidence lower limit on the ratio of up to down flux. As men- 
tioned above, due to geomagnetic effects this ratio is expected to be less than one due to a 
small enhancement of the downward flux. We neglect the effect of the geomagnetic enhance- 
ment of the downward flux in calculating our neutrino oscillation limits. This makes the 
limits a bit conservative since any reduction in the upward relative rate due to geomagnetic 
effects will instead be ascribed to possible oscillation effects. 

In the region 2.5 x 10~ 4 < Am 2 < 7.0 x 10~ 3 eV 2 where the upward neutrinos have 
traveled over several oscillations and the downward neutrinos have not had a chance to 
oscillate this bound on the ratio can be converted into a limit on sin 2 (29). 

1 - ^ sin 2 (20) > 0.64 
sin 2 (26>) < 0.72 

This is in direct conflict with the results || sin 2 (29) > 0.82 within the mass range 5 x 10 -4 < 
Am 2 < 6 x 10~ 3 eV 2 . A possible source of this discrepancy is a misinterpretation of their 
angular distribution in terms of low mass neutrino oscillations || by the Super Kamioka 
collaboration. 

A plot of the region excluded by this analysis and the analysis of reference [[| is shown 
in figure p]. The section marked "Super Kamioka Allowed" on the figure is not their true 
allowed region but an outline of the maximum range permitted for some values of Am 2 and 
some values of sin 2 (2$). The Super Kamioka fit implies that sin 2 (29) = 1. Figure [l] was 
calculated by integrating over the observed energy and distance distribution. 

This rate analysis fails to exclude Am 2 > 2.5 x 10~ 2 eV 2 because the upper hemisphere 
starts to show evidence of oscillations (in this low energy sample) so the downward rate 
would also be reduced. A shape test could be used to extend the results in this area. (The 
fact that the upward rate is not significantly greater than the downward rate permits one to 
exclude the range Am 2 < 0.1 eV 2 . To calculate detailed limits one must integrate over the 
source to detector distance since the muon decay length is comparable to the path length.) 
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The region 6.3 x 10~ 5 < Am 2 < 1.0 x 10~ 4 is more sensitive to rate effects than our simple 
analytical analysis above indicated. In this mass region the oscillation hypothesis takes the 
bulk of the upward data through the maximum. The ripples in the figure are indications of 
the repeat of this effect for the second and third oscillation etc. 

These results confirm the excluded region calculated from a study of an independent 
sample of upward going muons by 1MB JF|. 

A larger sample of 1MB 3 contained event data || is available and could be used to 
extend these results. With double the data sample the upper bound on sin 2 (2$) would drop 
to 0.54. The Super Kamioka single ring showering events H should yield a limit on v e — > v e 
of the order of sin 2 (2# ea .) < 0.23 over a comparable range of Am 2 . Precise contours depend 
on details of the measured v e energy distribution. 
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FIG. 1. 90% confidence level excluded region 
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